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The gravitational-wave instability of r-modes in rapidly rotating compact stars is believed to spin them down 
to angular frequencies Q. ~ 0.1flK ep ier soon after their birth in a Supernova. We point out that the r-mode 
perturbation also impacts the neutrino cooling and viscosity in hot compact stars via processes that restore weak 
equilibrium. We illustrate this fact with a simple model of spin-down due to gravitational wave emission in 
compact stars composed entirely of three-flavor degenerate quark matter (a strange quark star). Non-equilibrium 
neutrino cooling of this oscillating fluid matter is quantified. Our results imply that a consistent treatment of 
thermal and spin-frequency evolution of a young and hot compact star is a requisite in estimating the persistence 
of gravitational waves from such a source. 
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I. INTRODUCTION 

Neutron stars that undergo pulsations of quadrupole order or higher can emit gravitational waves. Such quasi-normal pulsation 
modes are typically expected in the aftermath of a supernova or a sudden rearrangement of the neutron star crust or core (egs., 
a starquarke or a hadron-quark phase transition). The possibility of detecting the associated gravitational radiation in current 
experiments such as LIGO and VIRGO is enhanced by the existence of the r-mode instability in rotating neutron stars. Discov- 
ered by Andersson |Q]], and detailed in subsequent works J2,[3t], r-modes are fluid pulsations whose dynamics are governed by 
the Coriolis force, and the I > 2 modes grow by emitting gravitational waves in all rotating neutron stars irrespective of their 
rotation speed. So long as viscous damping effects are small, which is the case for hot young rapidly rotating neutron stars, 
the r-mode can become unstable to the emission of gravitational waves while the star spins down by losing angular momentum 
(the Chandrasekhar-Friedman-Schutz or CFS instability HHH]). This is one of several reasons why neutron stars are regarded as 
promising targets of opportunity for gravitational wave detection |@]. 

For experimental detection of gravitational waves (see eg., ^ for a pedagogical treatment), it is important to have theoretical 
predictions for the gravitational waveform (eg. h{f) in the frequency-domain). According to the simple model discussed in J3] 
which we employ in this work, the angle-averaged signal in the time domain h{t) depends on the distance, instantaneous angular 
velocity £2(f ) and average density of the star, as well as the r-mode amplitude oc(f ). The time evolution of the mode amplitude 
and of the angular velocity constitute two coupled first-order differential equations which can be solved explicitly for a particular 
stellar configuration once viscous damping and gravitational wave timescales are known. Since the viscous damping timescale 
is very sensitive to the temperature of the fluid, thermal evolution of the star affects r-mode evolution. In previous works J3- 
IH, while this fact has been recognized, the converse effect, viz., that the r-mode perturbation also affects the cooling rate of 
the neutron star has been missed. A typical r-mode perturbation induces chemically equilibrating weak reactions that alter 
the neutrino emissivity from the standard equilibrium expressions and hence impact the cooling rate quantitatively. Therefore, 
adopting a simple power law for the cooling rate based on the temperature-dependence of equilibrium emissivities may not 
accurately describe the spin-down of the star. Thermal evolution, r-mode evolution and spin-down (i.e, rotational evolution) of 
a young and hot compact star are all coupled and a consistent analysis of these must be performed in any study of persistent 
gravitational waveforms from a rapidly rotating compact star. 

In this paper, we illustrate this generic interplay of thermochemical and spin-down effects and quantify it in a particular case, 
viz., that of r-modes in three-flavor (u,d,s) degenerate quark matter. Such a phase can arise in the ultra-dense interior of a 
neutron star or comprise almost the entire matter in a quark starQ. We also take into account leptonic (Urea) contributions to the 
viscosity of strange quark matter, calculated recently fill . Our main findings are: (i) the star can undergo two distinct epochs 
of gravitational wave emission and rapid spin-down, of very different duration, separated by a decades-long hiatus where its 
rotation rate is stable; (ii) the leptonic viscosity can be important for the spin-down history of the star if the initial spin-frequency 
V < 600 Hz (v w 0.5£2{Cepier/ (2ft)); (iii) the star takes about an order of magnitude shorter time to spin down to a stable frequency 
if non-equilibrium neutrino emissivities are considered and the evolution equations are evolved consistently. 

In Section [II] we provide a general overview of the link between r-modes, viscosity and gravitational waves, outlining the 
phenomenological model for mode and spin evolution introduced in J7[| and adopted for this work. In Section [III] we explain 
how r-mode perturbations affect the viscosity and neutrino emissivity of the constituent matter. In Section llVl we present results 
from a coupled analysis of thermal and spin-down evolution of a strange quark star. Section [V] gathers our conclusions and 
discusses the implications for persistent gravitational waves from neutron/quark stars. We include some mathematical details in 
the Appendices. 



II. fi-MODE EVOLUTION AND GRAVITATIONAL WAVES 

Pulsation modes of compact stars are classified by the nature of the restoring force in the perturbed Euler equation. Analysis 
in the co-rotating frame shows that the Coriolis force provides the restoring force for small non-radial fluid perturbations IU6I1 . 
It is conventional to assume that purely toroidal mode solutions to the perturbed Euler equation exist even for rotating stars, in 
which case the velocity perturbation can be expressed to (i2) as ^ 

hv=R£lf{r/R)YfJ^' , (1) 



Although cold and dense quark matter is believed to be in a color superconducting state fToh . we will not consider additional facets introduced by Cooper 
pairing here, leaving such considerations to future work. For r-modes in color superconducting quark matter, see HlHl4ll 
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with R and £2 being the star's radius and angular velocity, / and Yf m the radial and angular functional dependence of the 

mode, and the mode frequency in the co-rotating frame. If we restrict ourselves to isentropic perturbations (no composition 
or temperature gradients), it can be shown that only the I = m modes survive 11711 . that f(r/R) = a(rfR) 1 with a an arbitrary 
constant (dimensionless mode amplitude), and that lfl8ll 



^ = -^- + 0(tf). (2) 
(m+1) 

These are the r-modes and they are termed quasi-normal when they lose energy to gravitational waves (/ > 2). The energy of 
the r-mode is dissipated according to JH]: 



dE 

-d7 = - {& ° 



-m£2) 2 " , + 1 co rot |8/ mm | 2 -| £ / 3 r(2ii8o" i 8o;„ + C8o8o*) , (3) 



where 8a = V„8v a is the volume expansion due to the r-mode and 8o„/, = j(V a 5vt, + Vi,8v„ — |8„/,V c 8i' c ) is (the traceless 
part of) the shear tensor. The first term is the energy radiated in gravitational waves to lowest order in £2 with 8/ mm being the 
current multipole 



S/ mm -fl 2 £2 J drp ( -j . (4) 

For m > 2, a) rot < »z£2, so that the r-mode energy grows with gravitational wave emission, triggering the r-mode instability. 
Viscosity from a variety of mechanisms in the fluid can suppress r-mode growth, and counter the instability. The fate of the 
r-mode is determined by a competition between the timescale for viscous damping and gravitational wave emission. According 
to the simple model for mode growth in the linear regime described in 0], the evolution equations for a (mode amplitude) and 
£2 (angular velocity) are: 



a = - 



t v l+a 2 Q 
a a 1 - a 2 Q 
Igr iv 1 + a 2 Q 



where the viscous damping timescale is Tv = (t^ 1 ') 1 (^: bulk viscosity, T|: shear viscosity) and the gravitational wave 

growth timescale is Tqr. In general, Tj" = —(dEj/dt) /2E; i = C,,r\,GR. Reading off (dEj/dt) fromEqn. (0, and using the lowest 
order expression for the r-mode energy E, one obtains the expressions for the damping timescales x provided in Appendix B. 
The quantity Q = 3 J/21 where 0] 



1 r R Sir r R 

J=—J pr 6 dr; 1= / pr * dr , (6) 

MR 4 Jo H 3MR 2 Jo 

with M being the stellar gravitational mass and p the energy density. In this model 0], exponential mode growth is followed by 
non-linear saturation characterized by a 2 ~ 0Cq ~ 1 where d w 0. While such large saturation amplitudes are probably unrealistic 
due to multi-mode coupling and onset of non-linear effects SEIli, we prefer to use this large value to demonstrate the 
qualitative effect induced by the coupling of temperature and mode evolution. Ultimately, our analysis should be applied to 
more realistic saturation amplitudes and include the coupling to daughter modes HUcJ. Within our working model 01, once 
saturation is achieved, only the angular velocity evolves 



2Q. a 2 Q 
Xgk 1 - alQ 



^=—^T^- (7) 



It is important to note that Xy is strongly temperature dependent while Igr is temperature-independent so that solution of 
Eqn. (0 requires an estimate of the thermal evolution T{t), which is determined in the simplest case (ignoring thermal gradients 
and conduction) by 
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dT 

C v — = -L v +L hl (8) 
at 

where C v is the specific heat, L v is the neutrino luminosity and L/, is the heating rate due to dissipation^. While the r-mode is 
growing, the local baryonic number density fluctuates according to 



n B (t) = n° B + 8n B (t) = n° B + 8« B Re(e ! '™) , (9) 

where co = co^ and hn B (proportional to a, see Eqn. ( TTZl i) is the magnitude of the fluctuation. These fluctuations, which 
first appear at o(Q. 2 ) rather than 0(£2), destroy weak equilibrium, considerably enhancing the neutrino emissivity above its 
equilibrium value. They also affect the bulk viscosity, so that ^ and L v must be evaluated in the fluctuating foreground created 
by the r-mode. In the next section, we describe this computation, highlighting the effect of the r-mode perturbation on ^ and L v . 



III. tf-MODES, VISCOSITY & NON-EQUILIBRIUM NEUTRINO RATES 

The r-mode perturbation described by the density fluctuation Eqn. dT2T > causes a change in the fluid volume as well as the 
chemical potentials of the electron (e) and the up (m), down (d) and strange (s) quarks. Since weak equilibrating processes are 
relatively slow, the pressure fluctuation, expressed as 8P = 8Poe ,t0 ', lags the density fluctuation through the former's dependence 
on the fluctuating chemical potentials (i.e., 8Po acquires a phase). The dependence on temperature fluctuations is ignored since, 
as mentioned before, dissipated heat can be carried away much faster than weak equilibration timescales. Both leptonic and non- 
leptonic processes in dense quark matter can contribute to the viscosity at temperatures T > 10 10 K for typical r-mode frequencies 
CO ~ 1kHz [ 15]. Following the approach of lfl5ll . it is useful to choose the two independent fluctuations in the chemical potentials 
corresponding to the leptonic processes as 8/ji = 8(/j^ — /j„ — p e ) and 8/J2 = 8(/j v — /j„ —He) which can be expressed in terms of 
variations in the following quantities: baryon number (8ng /n B ), charge (bX e ) and strangeness (5X S ) 

dni=A i 8X s +B i 8X e + C i (—); X s = — , X e = — , (10) 
\n B J n B n B 

where n s and n e are the strange quark and electron number densities. The non-leptonic process d + ima + sis characterized 
by 8/J3 = 8/J1 — 8/j?. The bulk viscosity, which is a measure of the power dissipated per unit volume over an oscillation cycle is 
given by 

1 «° 

C = Im(8Po)^-; Im(SPo) = ~ nB [CiIm(8Xe,o) + (C 2 - Ci)Im(8X Si0 )] . (11) 

CO 0«B 



For the density perturbation, we adopt an expression derived in 12 111 (also Eqn. (6.8) of 122I0 

^l=R^[8U Q + m, (12) 
n B dP 

Wo = a P£ +i (cos 9)e im *; 84> = aS^ («« + i (cos 9)e""* , 

where 5<t>o is the change in gravitational potential due to the r-mode oscillation (80o(K) is obtained from Eq.(4.4) of l22Tl ). 
The perturbation, and hence 8/j,, has a radial and angular dependence, leading to different emissivities and viscosities in different 
directions and radial distances. We take a much simplified approach whereby, as is common in interpretation of neutron star 
thermal observations, the temperature T refers to or is related to a local surface temperature map (eg., a "hot spot") l23ll . Even 
if full angular dependence is retained, additional anisotropy in the bulk viscosity can be induced by strong magnetic fields in 
quark matter J24J. For now, we proceed with our simplifying assumptions of homogeneity and isotropy with magnetic fields set 
to zero. These assumptions should be relaxed in fully 3-dimensional calculations, which require extensive computations on a 



2 Viscous heating by shear viscosity, estimated as in 01, has a small effect on the temperature evolution in case of strange quark matter, while local heating due 
to bulk viscosity is rapidly radiated away in neutrinos and is unimportant for even the smallest timescales required for evolving eqs j5J. Therefore, we only 
include viscous heating from shear viscosity in . 
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spatial grid. The computation of Im(8X e)0 ), Im(8Xj o) an d the co-efficients Ai,Bi,Ci is straightforward once an equation of state 
is specified (see Appendix A for the equation of state for degenerate strange quark matter adopted in this work). We find 



^ = 7! 72^2 E Wflij, (13) 

W Si +82 (ij)=l;i<j 

gi = (AiB 2 -A 2 Bi)£-II; g 2 = aia2(A 2 -Ai) + aia3(A 2 -B2) -a 2 a 3 Bi , 
E = ai + 0C2 + 0C3; n = aia 2 0C3 , 

012 = £(AiC 2 - A 2 d f + n(Ci - C 2 ) 2 , 13 = E [Ci (A 2 - B 2 f + C 2 B 2 ] 2 + ncf , 23 = E(Ci2» 2 - C 2 fii ) 2 + flC^, (14) 



where A, ,B, ,C, are given in Appendix A. In Eqn. ( fT3l , a, = «b(o/A,, where the A.,s are obtained from the decay rates (F) of the 
non-equilibrium processes at first order as 



-^net = r c /->M+f-+v f ~^u+e-^d+v e ~ A-18/J1 , (15) 

^net = ^s->u+e~+\ e ~ ^u+e-^fS+Ve re ^2^2 , 

Tnet = Tn+^i+H _ r s+ „^„ +f / PS A,3(8^l - 8^ 2 ) = A-38^3 , 

where the second relation in each line is a commonly made approximation to first order in 8/7. The net decay rates for the 
leptonic processes can be computed exactly from the corresponding non-equilibrium neutrino rates f25ll 



^-3(rS?-lf)=3r&; i=l,2 (16) 



35p 

where T/ and Tf, are the forward and backward reaction rates respectively, and e[g t = el = 6^, + By.,. As shown in 

Appendix B, the relevant non-equilibrium neutrino rates are simply related to the equilibrium emissivities (Etot ) as lEi 



,(i)>eq 



10nuf + 315uf + 21uf 
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(17) 



Combining Eqns.dT5li.dT6l and dTTJ with Eqn. (|2TT > generates a self-consistent solution for the 8/j, even for large m,-. Since 
it turns out that 8/j, ^> £!T can occur for the r-mode perturbation, we retain the full expression in Eqn. ( fT71 ). Inclusion of these 
"large-amplitude" rates has been recently emphasized in the context of r- modes [ 27 ] . The decay rate for the non-leptonic process 

r^j has been calculated in various works. We adopt the analytic result of Madsen l28H29ll . which improves Sawyer's result [f30h 
by including a 8/J3 term as well as m s ^ corrections (their inclusion increases the rate). This approximation matches quite well 
with the exact numerical result up to T ~ 10 n K, which is the upper temperature bound considered in this work. The equilibrium 
neutrino emissivities appearing in Eqn. (T7\ are lf3~lTl 



457 

e tot' eq = T^a s T 6 G 2 F cos 2 QcMeMdMu , (18) 

(2),eq_ 45771 6 2 2 , 

e tot - TggTJ^ Cr F sin V c fi s m s . 



IV. NEUTRON STAR COOLING AND SPIN-DOWN 



As mentioned in the previous section, 8/j,- can be much larger than kT over some time domain of the evolution. In our 
simulation, this happens in the initial few hundreds of seconds of the star's spin down due to the r-mode instability. Consequently, 
non-equilibrium neutrino emissivities are much higher than equilibrium values, and we should anticipate a change in the cooling 
history of the star. In Fig Q] (bottom panel), we display results for the angular velocity versus temperature evolution from the 
combined analysis of Eqn.((5]l, for which we used the bulk and shear viscosity damping timescales displayed in the top panel of 
the same figure. Cooling with equilibrium and non-equilibrium neutrino emissivities (Eqn. dT7j) is studied. 
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Firstly, we note that the double-valley feature of the damping timescale (top panel) is reflected in the double-peak feature of 
the critical frequency curve, which is the locus of points satisfying l/tf(T)\o. c = [l/x^(T) + l/x^T) + l/icw(T))]n c =0. The 
critical frequency curve marks the boundary between rapid spin-down and unstable mode growth accompanied by gravitational 
waves in the region above the curve and stable rotation with a decaying r-mode amplitude below it. As expected, the critical 
frequency curve mirrors the behaviour of the bulk viscosity of strange quark matter (see eg., Fig. 3 of lfl5lo . The larger of the 
two peaks corresponds to the maximum in the bulk viscosity of the non-leptonic process 02811 while the smaller peak is from the 
maximum in the leptonic process lfl5ll . The peak height is inversely proportional to the r-mode frequency (co = 0.88kHz for this 
figure). For the simple model assumed here, these features determine the evolution of the angular velocity with temperature, as 
the star cools and spins down. 

To determine the persistence of gravitational waves from r-modes in a cooling compact star, we track the angular momentum 
evolution against the critical frequency curve to see when the star enters/exits the stable region. In the case of equilibrium 
neutrino rates, which ignores the effect of the r-mode on neutrino emissivities, we find that the star rapidly spins down in about 
a few hours (w 10 4 seconds) from its initial spin frequency, assumed to be the Kepler frequency (Qk/ (27t) « 1200 Hz) down to 
about a third of the starting value. This is the segment OA in Fig. Q] (bottom panel). The time to spin down to a stable rotation 
rate does not depend on the initial mode amplitude (typically, we varied this between 10~ 6 and 10~ 3 ), nor does it depend on 
the mode saturation amplitude oco, which we varied between 0.3 — 1.0. However, the stable frequency at which the star settles 
down upon entering the stable region depends on oco, becoming smaller with increasing oco. Once the star enters this stable 
region, where bulk viscosity is large enough to damp the r-mode, the mode amplitude decreases to negligible values. Since the 
r-mode is damped, it does not affect the neutrino rates in any way and equilibrium rates apply. While a small time step of about 
0.01 seconds is required in the initial rapid cooling phase, the time step in the stable region can be safely increased to hundreds 
of seconds, since temperature evolution and angular velocity evolution is extremely slow. The star takes of order 10 9 seconds 
(nearly 100 years) to cool down to the point where it can once again enter the instability region. This hiatus between spin-down 
epochs is depicted by the segment AA 1 . In the interim, the star has cooled by two orders of magnitude down to T w 0.7 keV, 
but hardly changed its spin-frequency (here, we do not take into account spin-down due to dipole braking or the possibility of 
glitches). Upon entering the unstable region at A', the mode amplitude begins a rapid growth from a negligibly small value and 
saturates within a year. At this point, the star spins down again due to the CFS instability and takes O (10 8 ) seconds (10 years) 
to reach a stable value of £1/ (2k) ~ 100 Hz. During this phase, the star is expected to once again radiate gravitational waves. 
Thus, a strange quark star might be expected to become unstable to r-modes over two distinct epochs, separated by about a 100 
years, with the second one lasting about 10,000 times longer than the first. 

For a consistent and more accurate estimate of the persistence of gravitational waves, non-equilibrium neutrino rates must be 
taken into account. A clear difference can be seen in the angular momentum history in Fig. J] Although we begin with the same 
initial conditions, neutrino cooling is more rapid, since the emissivities are larger (Eqn. (fTTIi). Consequently, the star enters the 
stable region at a higher angular velocity and in less time than the equilibrium case (segment OB). We find this time to be of 
order 10 3 seconds (less than 1 hour). The cooling in the stable region is once again marked by slow cooling and a negligible 
r-mode amplitude, resulting in equilibrium-type cooling. However, the fact that the star is at a higher stable rotation rate now 
implies that it spends less time in the stable region (segment BB'). This lasts approximately 5.10 s seconds (« 15 years), to be 
contrasted with the equilibrium case, where this era lasts for about 100 years. Upon exiting this region, the star spins down 
rapidly once more, taking « 10 7 seconds (4 months) to reach its final stable value. This is the segment B'O' , shown in enlarged 
inset view in Fig. [l] The final rotation rate is exactly the same irrespective of the cooling history. Therefore, the importance of 
incorporating non-equilibrium emissivities is only evident in the first 100 years or so of the star's spin-down history. Another 
important difference is the shorter time for which the star woud be active in gravitational waves if non-equilibrium rates are 
used. This indicates that the epoch of gravitational waves is shorter than expected, narrowing the window of detection. Similar 
results are expected for the case of neutron stars on the basis of the generality of the evolution equations, even though neutron 
matter has very different viscosities and cooling rates. 



V. CONCLUSIONS 

We have shown that the r-mode perturbation affects the neutrino cooling and viscosity in hot compact stars via processes that 
restore weak equilibrium. Using a simplified model of spin-down due to gravitational wave emission in compact stars composed 
entirely of three-flavor degenerate quark matter, we find that non-equilibrium neutrino rates associated to the r-mode fluctuation 
change considerably the cooling history of the star, as well as the angular velocity evolution, resulting in shorter epochs of 
gravitational wave emission for isolated quark stars. We have taken into account a recent calculation of leptonic contributions to 
the viscosity of quark matter, but find that it can only play a role if the initial spin frequency is chosen sufficiently small (« 600Hz, 
see Fig. IJi. This is because of the smaller viscosity of the leptonic process, specially at r-mode frequencies characteristic of very 
rapidly rotating stars. However, the role of the non-leptonic process is very prominent. It leads to two distinct rapid spin-down 
events in the star's evolution, caused by the CFS instability towards gravitational wave emission. Although we have focused 
on quark matter here as an exemplar, it should be clear that a similar exercise can be performed for any sufficiently simple 
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phase of dense matter, including beta-stable neutron-rich matter. It is important for r-mode saturation and damping studies that 
one includes the large-amplitude bulk viscosity results, as argued in |27|], as well as realistic saturation amplitudes [9]. While 
analytic approximations for the bulk viscosity exist for a simple quark matter equation of state like the one adopted in this work, 
a full numerical solution is required for hadronic matter, where non-linearities in the chemical equilibration rates can take more 
complex forms. Furthermore, non-equilibrium rates l32tl and updated viscosities l33ll need to be taken into account in performing 
a similar calculation with hadronic equations of state. Despite these complications, one can predict that for isolated compact 
stars made up of neutron-rich matter, since there is no stable region between « 10 7 — 10 10 K at rapid rotation rates J3, HI, 
non-equilibrium cooling is probably important only over the first year or more of cooling. 

We also mention some caveats that follow from the simplifications made in this work. The lowest order density fluctuation 
(Eqn. il2\ ) enters at order £2 2 , while the expressions for the r-mode frequency (Eqn.©) and the current multipole (Eqn. are 
(£2). Higher-order effects change the mode frequency by less than 10% (see eg., Table 1 in ifTHl ) while Eqn... suffices for 
gravitational radiation from the current multipole. A more detailed discussion of higher-order rotation effects is given in l22ll . 
Bulk viscosity timescales used in this work are approximate, and can be improved by replacing the Eulerian expression for the 
density perturbation, Eqn. ( fT2b by a Lagrangian one I22I1 . This can change the bulk viscosity damping timescale by an order of 
magnitude, but the critical frequency curve is hardly modified because the bulk viscosity has a steep temperature dependence. 
Thus, the results presented in the bottom panel of Fig.|T]are accurate in this respect, even though the bulk viscosity timescale, as 
shown in the top panel is only approximate. 

Non-equilibrium cooling can also be important for accreting neutron stars, which are believed to be persistent sources of 
gravitational waves (34j], undergoing spin-cycles as they get spun-up by accretion torques and spun-down repeatedly by the 
r-mode instability. Larger duty-cycles are expected for exotic phases, such as quark matter or hyperonic matter j35il, provi ding 
the motivation to extend our work to low-mass X-ray binaries. Even for isolated strange quark stars, Andersson has shown Dal 
that if a strange star enters the instability window spinning near the break-up limit, it will emit bursts of gravitational waves 
as the spin evolution straddles the critical frequency curve. We plan to explore this interesting case as well as the spin-cycles 
of LMXBs including non-equilibrium cooling rates and improved quark matter viscosites. Finally, we note that LIGO data 
sets ll37ll can already be adapted to search for imprints of the r-mode instability in gravitational waves 13811 . In anticipation of the 
improved sensitivity of advanced LIGO, it is important to continue to examine various factors affecting the r-mode evolution. 
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APPENDIX A 



With reference to Eqn. ( fT3l , the definitions 



ei =E(BiC 2 


-B2C1); 




fx = E(CiA 2 


-AiC 2 ); 
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g x =Z{A l B 2 


-B1A2)- 


n; 



lead to the expressions for hX s $ and 8X e .0 



(19) 



8X sfi = e -l±^(^); 8X e ,o = ^(^), (20) 



H + 'g2 V n B J gX + lg2 \ n B 

from which it follows that 



C/COS (cof ) + . 



-^A/cos (cof - - 

-si 



Iff- 



f- 

|-B,cos (o)t-§f) 



(21) 



with 
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tan^ = eig2 ~ e2gl ; tanfr, = /lg2 = /zgl . 
<?igi+e2,?2 flgl+hgl 

Given the following thermodynamic potential for 3-flavor degenerate quark matter (with electrons) 



(22) 



u, d, s.e 



12k 2 



Q 



4 

47t 2 

1 + 



i-4(i-|4)+^>ra] 



2a, 



(z}\n[F(z)} y / l _ 4) 2 -2(l -4) 2 J | 



-; Zf = mf / Mfj a s = g 2 /4n. 



where / denotes flavor, the coefficients Ai,Bi,Cj follow from 



(23) 



A 3 « O 

Ai = -n B Tr — ; B\ = -n B 
on d 

dfi s 

M=riB^—; B 2 = -n B 
upon using «, = — and the expressions 
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; c 2 = 



Ci = 



«J5 '^3 — 



3/*j _ 3/*« _ 3/4 
arc* dfZ(, an e 



(24) 



drif 3 \ ^ Jt f 



3/i e _ H e 
dn e 3 



(25) 



APPENDIX B 

Neutrino emissivity: Omitting volume normalizations which ultimately cancel, the emissivity is given by 



n4)4 (!,ftE "- ir|i ' |P ' 



spin 



(26) 



where (pi,Ei) denote the momenta and energy of the incoming and outgoing quarks and (co^, k) label the neutrino energy and 
momenta. The delta functions account for total energy and momentum conservation, and the function T (E ) = Y\in f{Ei ) Flout ( 1 — 
/(£,)) is the product of occupation (incoming) and blocking (outgoing) statistical factors for the quarks, with Fermi-Dirac 
distribution function /(£,) = (exp(p£' I ) + 1) (the neutrinos free stream). The matrix element rW is calculated from the 
corresponding tree level Feynman diagrams (see lfl5ll ) in the low-energy limit. It can be shown that for the non-equilibrium 
case (Ej w //, + 8/i,-), the neutrino energy integral in the total emissivity expression is modified from the equilibrium case l3lll 
such that 



m = T 



% 2 +x 2 



1 +e v 



kT J \ IcT 
yielding Eqn. (|17t . 

Damping timescales: For viscous damping of the r-mode, the shear viscosity damping timescale is 



(27) 



1 _(»-l)(2» + l) ,* ^ 



J,ft/rpr 2m + 2 JO 

where the shear viscosity is taken from l39ll . while the bulk viscosity damping timescale follows from 



(28) 
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1 1 fdE\ 

~~ ~2E \dj) 



(29) 



2m 

The gravitational radiation growth time scale is J2l 



/•« / r \ 2m+2 

E = ^ 7 _{m+iy{2m+l)iR 4 Q. z J drp 



1 327cGn 2 '"+ 2 (/n-l) 2m ^m + 2\ 2m+2 ^ J 



Tew c 2 ' n+J [(2m+l)!!f V m + 1 



/ ofrp 
70 



r /m+z . (30) 
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FIG. 1. Top panel: The temperature dependence of the shear (r)) and bulk (Q viscosity damping timescales from Eqs. ( 128b . J29l >. The stellar 
gravitational mass and radius (Newtonian gravity) is 1.49M® and 9.47 km respectively, while the r-mode frequency is 0.88 kHz from Eqn. l(2j- 
Bottom Panel: Spin down history versus temperature for equilibrium and non-equilibrium neutrino cooling (see text in Section llVl for details). 



